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Abstract 

We analyse de fluctuations of charge of the (1+1)D Dirac's fermion with charge conjugation breaking. This 
is done taking the separation between background soliton and anti-soliton going to infinity. 

^: 

■ 1 Introduction 

The interest in the investigation of non-integer charge has become stronger in quantum field theory since the 
04 ' original work [I]. This happened because this charge conjugation symmetric (l-l-l)D model, with a soliton anti- 
soliton background, predicts | as an eigenvalue of the charge operator and has been used to describe properties 
of the polyacetilen[21 The concept of fractionary charge has been reinforced when it was shown that there is 
, no charge fluctuations [H [5]. These ideas opened a door to investigations in condensed matter using, as a tool, 
' quantum field theory. 

The possibility that some system could have non-integer charge was also investigated in a variety of other 
contexts. Excitations of charge condensates were studied by Rice at all More recently, there has been renewed 
interest for non-integer quantum numbers. The construction of high quality samples of grapheme made it possible the 
discussion of two-spatial dimensional field theory models applied to these materials. Distortions of the honeycomb 
lattice of grapheme lead, in the continuum, to theories with fermions coupled to scalar fields (complex Higgs fields) 
and their topological excitations such as vortices induces mid-gap states. This state of affair reproduces much of 
the excitement previously experienced with polyacetilene in this higher dimensional analog. In particular, charge 



' fractionalization through the breaking of energy-reflection symmetry reappears [2l IHl E] • This mechanism should 
be contrasted with the previously studied phenomenon of fractional charges as in the fractional Hall effect, in two 
spatial dimensions , which rely on the forced breaking of time-reversal symmetry associated to topological order. 
This new mechanism is much closer to the ones witnessed with polyacetilene. Instead of fractional charges generic 
values for the charge are obtained. 

It is very natural that variations of the original models predict a much richer spectrum of possibilities. The 
simplest generalization is to include breaking of charge conjugation in the Jackiw's first proposal [3] as done in [TU] 
who obtained general non-integer charges. Further investigations have been done [111112) . In [T3] a general method 
for computing the charge has been applied to two- and four- dimensional models. 

However it is missing, up to now, the investigation of the fluctuations of the charge within models without charge 
conjugations symmetry. To perform this investigation within the simplest case is the goal of this work. 

In section 2, we present the model and construct the localized charge operator. In following section we compute 
the charge eigenvalue. The fluctuations are computed in section 4. 
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2 The number operator 

Let us start with the (1+1) dimensional Dirac equation 

£;^' = {(T2P + ai(f> + 7(73)^'. 
In terms of components, the equations of motion are 

{-ip + (t))ip2 = {E~j)^pi, 
(ip + 0)Vi = {E + j)^P2. (2.1) 

Here p = —id/dx, 4'{x) presents an soliton-anti-sohton background profile 

( +00 , a; < 
(j){x) -00 , < a; < L 

[ +00 ,L <x 

and 7 is a constant parameter responsible for the breaking of the Hamihonian conjugation symmetry. 
The continuum energies are given by 



and the positive energy solutions can be written as 



d+uk{x) 

d~Vkix), 



while the negative ones as 

^-(x),-( ^""^-(^^ ^ 



where 



--{k) 



\E\±l 



E 



It is remarkable that there is no zero-energy continuum, i.e., \E\ > + 0o- -f^'" the bound states, the analogous 
occur, i.e., the positive energy solution is given by 

' \ 9gVB{x) 



and the negative one is written as 



9gUB{x) 

'0%vb{x) 



where the parameters 0^ are defined in terms of the the bound-state energy in analogy to definitions of 0^(k). 
Note that thes parameters are equal to 1 when 7 is iqual to zero, so that the functions u's and v's become the 
upper and lower compontents of ^ . In this presentation the functions Uk{x), Vk{x), ub{x) and vb(x) satisfy the eq. 
(|2.ip with 7 = 0. All contribution of the conjugation breaking is codified by the parameters 9"^. The bound states 
energies are Eb = ±1/— + 0q + 7^, where the parameter n is defined by a transcendental equation which does 
not depend on 7, 

(jy^ — — 00 exp {~2kL). 

The normalizations are 



2 



/ uluk= / 

Jx J X 



Uf.Uk 

UgUB 



VkVk 



-Mk-k% 



VbVb = 2 



(2.2) 
(2.3) 



and the overlaps between ub and Uk as well as between vb and Vk vanish. 
Observe that the "conjugation" transformation 



* — > (cias + C2)*^ 



with 



and 



1 / 61+^61-^' 



e+9- 



1 



9+8- 



maps positive to negative energy states. It is important to note that, in spite of the conjugation mapping 
7^ Pk"^^- present case, defining the first quantized scalar and chiral densities p = and p = 



(2.4) 
(2.5) 

(2.6) 

; above. 



~E>0 



ptt'' = {ci'+c,^)pr^>+2c,c,ru 

This implies that the fermion charges acquires a generic value, not restricted to half-integer values [?]. 
The second quantized quantum field operator can be expanded as 



-iEkt 



0+Ukix) 



ak + e 



e-ui{x) 
-e+vi{x) 



with canonical anticommutation relations between a^, 6fc, gb and 6b, i.e. 



{«fc.4'} 

If we define the localized number operator as 



[hk,hl,] = 5{k-k'), 



Nf = 



where f{x) plays the role of an smearing function, peaked near the soliton, we obtain 



(2.7) 



(2.8) 



(2.9) 



Nf 



x,k,k' 



which can be expressed as 



ilak{e+e'+ui,uk + e e' vi,vk) - blbk{e e'~ui,uk + ^ 



+e-^(E.+E„)t 



Uk 



0-0'+ Vk' 



Vk'Vk) 



+e 



+i{Eu+Ey)t 



albl{0-e'+uk'*ui - e+e'-vk'*vi) 



+ls{k - k') [{6-^ - e+^)uluk - {0-^ - 0+^)vlvk 



Nf = Ni + Nl + Nl. 



(2.10) 



(2.11) 
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Here represents the quantum contributions of the bound state sector, Nj^ the continuum sector and Nc is the 
c-numbcr contribution coming from both the bound as well as the continuum sector after subtractions of the zero 
soliton background. 



Ni = 



J X 



Ofvl) 



(2.12) 



/ 



+ /J[(^B'-^f)(«l-4) 



(2.13) 
(2.14) 



where the subscript 0, in 6^, means that the energies are the ones with no soliton background, but with 7^0 and 



H = K 



■N, 



f++ 



Here it is useful to decompose the continuum contribution in terms of the creation and annihilation operators 



N, 



and 



/+- 



Jx,k,k' L ^ 



'vl,vk)-h\,hk{0 9' ul,Uk 



vl,Vk) 



H-- = I 

J i£ jk jk 



(2.15) 
(2.16) 

(2.17) 



In the last equations, 6' = 6{k'). 

It is important to note that there is a contribution to the c-number term coming from the continuum states which 
appears when 7 7^ 0. It is not eliminated by subtracting the vacuum contributions computed in the zero soliton 

sector. Indeed, although wc have cxplicitcd these zero soliton regularizing terms they turn out to be identically 
zero, since = l^'^p. These regularizing terms are necessary if we compute separately the contributions from 
Uk and Vk- 

It is convenient to rewrite these equations by defining the symmetrized and anti-symmetrized combinations of 
the factors defined below 

= {e-{k)e+{k') + e+{k)e-{k'))/2, 



and 



With these definitions we obtain 



e'^^' = {9-{k)9+{k') - e+{k)9-{k'))/2, 

9f = {e+(k)e+{k') + 9-{k)0-{k'))/2, 
9^1^ = {9+{k)9+{k') - 9-{k)9-{k'))/2. 



1 ' 


' ^-i{Ek-E^,)t 




i x,k,k' 







(4,afe - bl,bk)0's \ul,Uk + vl,Vk) + (4'"fe + bk'bk)9^A i'^k''^k - vl,Vk) , (2.18) 



Ni 



x.k.k' 



(2.19) 
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and 

nI++^ I f e+^^^^+^^'Hl,al 9i'\ul,ul-vt,vl)~0^2\<'<+vl>vl) . (2.20) 



' x,k,k' 

On the next section we will investigate the fluctuations of the eigenvalues of this operator. 

3 Charge fluctuations 

Let us first remark here that as / — > 1, Nf becomes diagonal, with integer eigenvalues. This corresponds to the 
number operator of the model with the soliton anti-soliton distance fixed. In this case the orthonormality of the 
functions can be invoked, leading to 

Nl-^Nk= I {alau-hlhk), 

Jk 

N{ — >Nb= - b^B^B 

and 



Now let us turn to the study of the limit L — > oo followed by / — > 1. Since the Uk{x) and Vk{x) functions are 
independent of 7, we have 

hm / f{ul+v%) = l (3.21) 



and 



hm / /(4-4) = -i, (3.22) 



where L is the soliton-antisoliton distance. 
In the following we use the point-wise limit 



hm \u,\' = (3.23) 
hm \v,\' = \vi:r, (3.24) 

I lim Ufc|2 ^ hm (3.25) 

where u]^^ and vj:^ corresponds to the solutions, with 7 = and with only one anti-soliton and no soliton. See the 
appendix A for the proof of equations (|3.23l - [?.25p . 

The limit of the bound-state energies, in the long distance case, is characterized by 



and 



with 



lim Eb = 7. 



In this limit we obtain that 



and that 



lim 

/-)-i,L->oo y Vb 



ub \ I 

,1-S 



lim 9+ = V2, hm = 0. 



From this we conclude that there is one state with negative energy E = —7 in the limiting case. 

The situation where L ^ 00 followed by / — > 1 is caracterized by some further usefuU properties, i.e. 



5 



lim / f {ul,Uk - vl,Vk) = 0, 
lim / / {uk'Uk - Vk'Vk) = 0, 
lim / / (u*,Ufe + vl,Vk) = T5{k - k') + R5{k + k') 



and 



lim / / {uk'Uk + Vk'Vk) = T'S{k + k') + R'5{k - k'). 
Using these relations and the explicit form of 6^ we obtain 

(T{k){alak - blbk) + R{k)ial,ak - bl.bk)) , 



/+- 

k 



N, 



/,±± 



and 
so that 



As has been shown in the appendix B, we obtain that 



iV/^zitan-i(^). 



(3.26) 
(3.27) 
(3.28) 

(3.29) 



(3.30) 



7 

Note that the limit when 7 goes to zero leads to the — ^ value expected. The same result is obtained by considering 
the case when we start with only one anti-soliton background, see appendix B. 

It is interestingly to note that the investigation of the case 7 7^ easier to deal than the case 7 = 0. This occurs 
since there is no mixing between the creation-destruction operators for the Nb term in the limit L going to infinite 
and / going to one. In contrast to 7 = case, no Bogoliubov transformation is necessary. Let us compute the 
fluctuations to the state 



IV >=6t|o,0 >= |0,1 > 



Using eq. (??), we obtain 



lim < Vl-^V-'^IV' > 

/— >1,L— >-oo 



lim <tp\Nl\tp>+ lim < Vl-^V^IV' > + 1™ < tp\Nl,\tp > 



= + 1 + Ng". 
Let 's consider the quadratic number operator 

(Nff = {Nlf + {Nir + {N[,f + + 2N[,{Nl + iV^). 

So we have to compute 



(3.31) 



(3.32) 
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lim <i^\{Nlf\^> = lim <ij\{Nl'+- +nI'++ + Nl'-~fM> 



= lim <^|(7V,{--.iV,^'+^)|^> 

+ {0^^^f{ulul, +vlvl,){uuu^, +vuv^,)}. (3.33) 

The first term in the last equation is discussed in [4^ and is seen to be in this limit. For the two middle terms the 
integrals of the factors with negative signs between u*u and v*v vanish for analogous reasons. That the last term 
also vanishes comes from p.26ll3.29| which makes 6'^^'' — > 0. The result is 



Furthermore 



and 



Also 



lim < ) IV' >= 0. 



lim <V|(iV+)2|V;>=l 



hm {N[,f = {N^f. 

j^l.L— i-oo 



hm <i,\{Nm)\i:>=Q 

/— >-l,L— ^oo 



and 

hm <^|2iV^(iV/ + iV^)|V'>=27V^. 

Collecting all terms we obtain 

< i}\{N+ f\i! {N^ f + + 1 = {N^ + 1)2, 

so that 

< ^\5{N+f\'4> >= 0, 

which means that there is no fluctuations. 

4 Conclusions and Discussion 

In this paper, we have studied the charge eigenvalues and their fluctuations for a model without charge symmetry 
that describes fermions, in (1+1)D in a soliton anti-soliton background. This has been done in the context of second 
quantization extending the analysis of p2] to the case of soliton and anti-soliton. Although this model had been 
throughously investigated in the literature it was remaining to be shown that the variance of charge expected value 
vanishes, which is a necessary condition to imply that this charge is a good quantum number. We have obtained 
zero fluctuations as the result of our analysis. 

It is interesting to note that in the case 7 7^ there appears a C-number contribution to the charge operator that 
misses in the case 7 = 0. This contribution, which comes from the continuum states, turns out to be responsible for 
the existence of the continuum spectrum of the charge when 7 7^ 0, varying between ±i. If we take a look in the 
resultant charge, eq. (|B.ip . we identify a part that is i and another that is continuous. The second part collects 
the contributions from the charge conjugation breaking term. Curiously, in a certain sense, the treatment is easier 
when 7 7^ 0. This occurs since in the presence of a soliton and an anti-soliton there are two finite-norm fcrmionic 
states centered in the soliton or the anti-soliton positions. Their energies go to zero when the distance grows for 



7 



the case 7 = but remain distinct in the case 7 7^ 0. This makes the overlapping of the solutions in the infinite 
distance limit finite when 7 = and leads to the necessity of a Bogoliubov transformation that is not necessary 
when 7 7^ 0. 

The strategy of this work has been to inquire about the physical meaning of the charge induced by the soliton 
through the imersion of the physical situation in a broader one where a soliton and anti-soliton are present. Then 
the simplest case should be re-obtained in the limit when the soliton to anti-soliton distance grows infinitely. If the 
fluctuations of the charge operator defined around the soliton were kept finite this would be a strong evidence of the 
partial character of the non-integer charge. In this analysis the point-wise limits eqs. p.23l - 13.25|) have played an 
important role. They contain the information that those systems, soliton anti-soliton and only a soliton background, 
are equivalent after we take the limit when the distance between soliton and anti-soliton goes to infinity. We have 
obtained these relations in the sense of distribution theory and have taken strict control on the orthonormality 
character of the functions involved in their definitions along the limiting process. 

Lets us stress that the lessons of our analysis should be applicable to field theory models of graphene which 
present fermion fractionalization in 2-l-lD also without conjugation symmetry. 
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Appendix A 

In order to evaluate the C-number contribution to the number operator, we start introducing 



into (|2.12p . which leads to 



TV/ 



/ 

x.k 2 



27 



\m\ 



27 , 2 



(A.2) 



We are interested here in the limit L — > 00 followed by / — > 1. For the bound-state case, the function ub goes 
to for any finite argument. The only contribution comes from vb resulting in i in the limit of the last integral. 

A much more careful treatment is required to the continuum. Let us consider here the square profile functions 
for the (j){x) soliton anti-soliton configuration 



h ,0<x<L 



We can proceed as in [4l and construct one explicit solution Uk.L,<i>oi^)j with the requirements 



Wfc,L,0o(2;)wfc',L,0o(a;)da; = ^S{k - k ), 
1 



i>o — >o 2-^7r 



Akx 



and 



lim Ufc,L,0„(a;) = —= 

Lq ^0 l^TT 



Akx 



Indeed, for instance in the case a; < 0, we obtain: 



(A.3) 
(A.4) 

(A.5) 



Uk,L,4,o{x) Vke''''' +^ke 



—ikx 



(A.6) 



with 



1 . 



ikL 



ikL ^ 



(A.7) 



Here 



Vk 



1 



2V^' 

i 



1-iX 



2v^ [1 + 2A2 + 2A\/rTA2 cos(A:0] VTTA^ 



expi—ikL) H , 



exp{ikL) + 2A]. 



(A., 



Similar expressions hold for < x < L and for L < x. 

We are now able to analyse the limit L — > 00 taking into account that rapid oscilating functions involved can 
be seen as distributions theories. For instance, 
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L — >oo a + bcos{kL) _|_ ^2 

1™ "Z^^^'nl = T- " ■ (A.9) 
L^oo a + &cos(fcL) 5 6Va^ - 

After some computations, in this approach, we obtain 

lim |ufc,,.,^j2=^K-p, (A.IO) 

L >oo 

where uj^^ corresponds to the solution of the problem with 1-antisoliton profile, 

^ ' \ -00 ,a; > 

For the Wfc,L,0o(a;), similar results can be obtained. If we follow the same steps, we can see that 

lim |wfc,L,0(,(2;)P 7^ I lim Uk,L,ci>oix)\'^ ■ (A. 11) 

L — ^00 L — >oo 

Indeed lim^ >ocUk,L,it>oi^) results in solutions of the 1-soliton profile that are not orthonormalized. That limit of 

the square is not the square of the limit, while not entirely clear on physical grounds, is understandable since the 
product of distributions is not in general well defined. 

Let us face now the problem of solving (|A.2p . The convergence factor f{x) allows one to interchange the space 
integrals with the limit L — > 00. We obtain that 



lim 



N, 



f _ 



lim 
lim 



x.k 



lim 

L — >oo 



2^ 

\E,\ 



27 



\Ek\ 



\Ek\ 



27 /„*1-S„,1-S 



(A.12) 
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Appendix B 

If we remember that, to the case 7 7^ and a one-anti-sohton background profile, there is only the negative 
energy {Eb = —7), the bound-state field is 







-\/2vb{x 



The number operator becomes 



Nf 



t 

J x.k.k' 



[alak - blbk) ~ b^B^B^-\ 



e'+ 



fe.: 



\m\ 



(B.l) 



Now the c-number contribution can be computed for instance by using the explicit form of the continuum functions 



,i-S|2 



(B.2) 



Using equations ()B.2[ IB.1|) and performing first the integration over x and recognizing that the oscilating terms 
in infinity do not contribute and then integrating in k we find that 



Nf = j{alak - blbu) - b\bB + i tan-i(^). 
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